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ABSTRACT 


This  report  is  concerned  with  the  dynamic  behavior  of  a  rocket 
motor  and  thrust  stand  combination  when  restrained  by  a  hydraulic 
force  balance  system.  The  system  is  analyzed  taking  into  account  the 
nonlinearities  caused  by  limiting  flow  and  fluid  compressibility.  The 
analysis  is  followed  by  a  discussion  of  the  data  obtained  from  an  analog 
simulation  of  the  system.  Investigations  of  the  basic  stand  without  the 
servo  are  followed  by  investigations  with  the  servosystem  installed  for 
comparison  purposes.  Both  the  time  and  frequency  domains  are  in¬ 
vestigated  using  ramp  and  sinusoidal  forcing  functions  of  thrust.  The 
results  of  this  investigation  show  that  many  desirable  advantages  may 
be  realized  using  hydraulic  force  feedback.  The  control  signals  may 
be  adjusted  to  provide  zero  displacement  of  the  rocket  motor  position 
during  the  steady  portion  of  the  firing  and  at  the  same  time,  provide  a 
high  degree  of  damping.  Concurrent  with  these  conditions,  the  accel¬ 
eration  of  the  rocket  motor  may  be  greatly  reduced  during  the  firing 
transient. 
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SECTION  I 
INTRODUCTION 


The  testing  of  turbojet  and  ramjet  engines  under  simulated  environ¬ 
mental  conditions  paved  the  way  for  similar  testing  procedures  with 
rocket  motors.  At  the  Arnold  Engineering  Development  Center  (AEDC), 
the  transition  has  involved  major  modifications  in  test  cells  along  with 
construction  of  new  ones.  Instrumentation  needs,  data  gathering  tech¬ 
niques,  data  reduction,  and  general  operating  procedures  have  been 
modified  and  improved  to  satisfy  the  requirements  of  modern  testing. 

From  the  standpoint  of  instrumentation  interests,  the  switch  to 
rocket  motor  testing  can  best  be  described  as  a  shift  in  emphasis  from 
steady-state  to  transient  phenomena.  Jet  engines  were  generally  oper¬ 
ated  at  fixed  environmental  conditions  for  each  data  point  and  ample 
time  was  allowed  for  all  transients  to  die  out  after  reaching  the  desired 
condition.  As  a  result,  most  of  the  transducers  and  recording  devices 
were  only  required  to  have  good  static  characteristics.  The  relatively 
few  special  tests  involving  throttle  bursts,  combustion  instabilities, 
ignition  and  flameout  transients  were  usually  of  low  frequency  content. 
Although  special  recorders  were  frequently  used  to  record  these  tran¬ 
sients',  no  stringent  dynamic  requirements  were  generally  imposed  on 
electrical  measuring  devices. 

The  criteria  governing  the  design  of  thrust  stands  for  jet  engines 
were  relatively  simple.  First,  the  stand  was  required  to  hold  the 
engine  in  place.  Secondly,  the  stand  and  load-cell  combination  was 
designed  to  have  as  high  a  natural  frequency  as  possible.  It  was 
assumed  that  low  frequencies  (in  the  flat  region  of  the  stand  frequency 
response)  would  pass  through  the  stand  to  the  load  cell  with  good 
fidelity.  Perhaps  the  most  serious  drawback  to  this  approach  was  that 
systems  of  low  damping  resulted.  Noise  could  easily  excite  the  thrust 
stand,  and  this  vibration  became  the  disturbance  source  for  instruments 
mounted  on  the  engine. 

Thrust  measurements  for  rocket  motors  placed  new  emphasis  on 
thrust  stand  dynamics.  Fast  thrust  rise  and  tailoff  transients  so 
greatly  excited  thrust  stands  that  little  information  concerning  true 
thrust  during  these  transients  could  be  determined.  Early  analysis 
techniques  involved  everything  from  hand  smoothing  thrust  recordings 
to  the  application  of  band  rejection  filters.  The  futility  of  these  methods 
quickly  accentuated  the  need  for  more  sophisticated  approaches. 
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The  near-linear  second-order  response  of  most  thrust  stands  pro¬ 
vided  a  basis  for  study  by  early  investigators  (Refs.  1,  2,  and  3). 
Several  methods  under  the  general  heading  of  “thrust  compensation1' 
were  introduced.  These  techniques  were  mainly  concerned  with  deter¬ 
mining  the  forcing  function  (true  axial  thrust)  of  a  known  differential 
equation.  Nonlinearities  such  as  changes  in  mass  caused  by  burning 
and  nonlinear  spring  characteristics  could  be  included  provided  they 
could  be  defined.  Operational  amplifiers  were  used  for  computation 
and/or  simulation  and  each  method  could  be  used  on-line  if  the  neces¬ 
sary  preparations  were  made. 

The  thrust  compensation  schemes  referred  to  above  were  concerned 
with  dynamic  compensation  for  axial  thrust  measurements.  The  prob¬ 
lem  of  multicomponent  thrust  measurements  is  much  more  complicated 
since  thrust  measurements  (static  and  dynamic)  in  one  direction  include 
components  of  thrust  in  other  directions.  If  only  static  measurements 
of  thrust  were  of  interest  these  interactions  could  be  decoupled 
mechanically  or  by  analytical  computation  schemes  provided  the  stand 
settled  down  after  the  initial  thrust  rise  (Ref.  4). 

The  axial  thrust  component  is  generally  much  larger  than  other 
components.  As  a  result,  measurements  in  the  vertical  and  side  direc¬ 
tions  will  contain  a  higher  degree  of  distortion.  Also,  dynamic  excita¬ 
tion  and  ringing  in  the  vertical  and  horizontal  planes  will  generally 
originate  from  the  axial  thrust.  For  these  reasons  a  high  degree  of 
rigidity  and  near-critical  damping  in  the  axial  direction  is  to  be  desired 
in  order  for  static  resolution  schemes  to  have  meaning. 

The  use  of  electrodynamic  actuators  as  active  force  feedback  ele¬ 
ments  to  produce  artificial  damping  has  been  investigated  (Refs.  5 
and  6).  Their  excellent  frequency  response  characteristics  when  used 
with  special  driver  amplifiers  make  them  well  suited  for  this  applica¬ 
tion.  As  force  actuators  in  test  cell  work,  however,  they  are  not  well 
suited  because  of  the  force  limitations,  required  cooling,  and  large 
physical  size -to -force  ratios. 

This  report  has  been  prepared  to  determine  the  performance  that 
might  be  expected  using  a  hydraulic  servosystem  to  alter  the  dynamics 
of  a  typical  thrust  stand.  The  study  involves  an  analysis  leading  to  a 
mathematical  model  of  a  thrust  system  and  an  analog  simulation  to 
study  its  behavior  under  controlled  conditions.  So  that  the  investigation 
will  have  practical  meaning,  typical  mass,  spring,  and  damping  coeffi¬ 
cients  as  experienced  with  actual  thrust  systems  have  been  used.  The 
specifications  governing  the  servovalve  and  piston  actuator  are  typical 
of  some  of  the  best  commercially  available  units. 
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Three  configurations  of  the  thrust  stand  and  control  circuitry  are 
considered  in  the  simulation  study  on  performance.  Simulation  data 
were  collected  and  are  discussed  concerning  both  the  time  and  frequency 
domains  for  each  configuration.  Configuration  1  deals  with  the  basic 
thrust  stand  without  the  servo  installed.  This  system  is  considered 
first  so  that  its  data  can  serve  as  a  reference  for  comparing  the  data 
taken  for  the  other  configurations.  Configuration  2  considers  the  same 
thrust  stand  with  the  hydraulic  servosystem  installed.  The  actuator  is 
mounted  on  the  thrust  butt,  and  the  piston  pushes  against  the  load  cell. 
The  servosystem  is  operated  as  a  positional  servosystem  for  this  con¬ 
figuration  to  position  the  piston  with  respect  to  the  actuator  housing. 

The  command  signal  is  furnished  by  a  manual  set-point  potentiometer. 
Configuration  3  is  identical  to  Configuration  2  except  that  the  command 
signal  is  furnished  by  live  measurements  from  the  load  cell  and  rocket 
motor. 

The  first  two  configurations  are  mainly  concerned  with  defining  the 
characteristics  of  the  thrust  stand  and  the  effect  of  adding  the  positional 
servosystem.  Configuration  3  is  mainly  concerned  with  the  composition 
and  adjustment  of  the  automatic  feedback  signal  to  provide  the  following 
characteristics:  high  degree  of  damping,  absolute  return  of  the  motor 
position  to  zero  after  the  firing  transient,  low  offset  of  the  motor  posi¬ 
tion  during  the  firing  transient,  and  a  large  reduction  of  motor  accelera¬ 
tion.  Visual  optimization  of  these  characteristics  is  accomplished  using 
the  repetitive  operation  feature  of  a  PACE  231R  analog  computer.  Data 
are  presented,  compared,  and  discussed. 


SECTION  II 

GENERAL  TECHNICAL  CONSIDERATIONS 


Before  a  detailed  analysis  of  a  particular  force  feedback  system 
can  be  undertaken,  it  is  important  that  the  stand  configuration  and  de¬ 
sired  mode  of  operation  be  understood.  So  that  this  investigation  will 
produce  qualitative  conclusions  concerning  the  advantages  of  the  method, 
the  stand  configuration  and  its  dynamic  characteristics  should  be  typical 
of  those  now  in  use.  Similar  practical  consideration  should  also  be 
given  to  the  hydraulic  servomechanism  and  the  measured  system 
variables  used  to  control  it. 

This  section  is  concerned  with  a  description  of  the  selected  thrust 
stand,  its  lumped  parameter  representation,  and  the  desired  perform¬ 
ance  under  controlled  conditions.  A  brief  review  of  the  general  fre¬ 
quency  response  characteristics  of  hydraulic  shakers  is  also  presented. 
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2.1  THRUST  STAND 

For  this  investigation  a  horizontal  engine  mount  and  thrust  stand 
system  is  considered.  The  configuration  is  shown  in  Fig.  la,  and  its 
diagrammatic  equivalent  is  shown  in  Fig.  lb. 


Piston  Load 

Actuator  Cell  Rocket  Motor 


a.  Schematic  Diagram 


b.  Equivalent  Circuit 


Fig.  1  Racket  Thrust  Stand  with  Hydraulic  Servo 

The  engine  is  represented  by  a  single  mass,  Mm,  and  the  load  cell 
by  a  spring -damper  combination,  Km  and  Cm.  Although  the  thrust  butt 
is  usually  a  complex  structure,  its  response  can  be  closely  approxi¬ 
mated  by  a  second-order  system.  For  this  reason  the  thrust  butt  is 
represented  by  the  spring,  mass  and  damper,  Kb,  Mb,  and  Cb,  respec¬ 
tively.  The  equivalent  mass,  Mb,  will  also  be  assumed  to  include  the 
mass  of  the  actuator  housing  and  servovalve.  The  mass  of  the  piston, 
however,  will  be  considered  separately  in  the  analysis.  The  displace¬ 
ments,  Xm,  Xb.  and  Xa,  are  measurements  with  respect  to  prefire 
locations  of  the  motor,  piston,  and  actuator  housing.  So  that  qualita¬ 
tive  comparisons  of  data  can  be  made,  the  thrust,  T,  will  be  a  limited 
ramp  function.  The  rise  time  and  maximum  level  will  be  varied. 
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2.2  PERFORMANCE  GOAL 

The  primary  purpose  of  the  hydraulic  servomechanism  is  to  pro¬ 
duce  the  necessary  expansion,  Xb-Xa,  in  the  thrust  link  to  hold  the 
rocket  motor  at  its  original  prefire  position  (Xm  =  0).  Also,  of  major 
importance  is  (1)  the  peak  value  of  Xm  during  the  firing  transient, 

(2)  time  required  for  Xm  to  reach  zero,  and  (3)  the  manner  of  return¬ 
ing  to  zero  (damping). 

On  the  basis  of  the  above  requirements  the  control  signal  to  the 
servoamplifier  consists  of  varying  amounts  of  the  following  signals: 

(1)  load-cell  output  signal  for  speeding  up  the  recovery  time,  (2)  integral 
of  Xm  to  make  Xm  in  steady-state  identically  equal  to  zero,  and 

(3)  velocity  and  acceleration  for  damping  control.  The  coefficients  for 
each  of  these  signals  are  determined  by  optimization  using  the  repeti¬ 
tive  operation  feature  of  the  computer. 


2.3  HYDRAULIC  SERVOSYSTEM 

Hydraulic  shakers  are  specially  designed  servosystems  for  pro¬ 
ducing  artificial  vibration  environments.  They  are  ruggedly  constructed 
and  capable  of  producing  large  forces  at  relatively  high  frequencies. 
Because  of  these  characteristics  the  specifications  of  a  commercially 
available  hydraulic  shaker  were  used  in  this  investigation  in  order  to 
arrive  at  realistic  values  for  the  servosystem. 

2.3.1  System  Description 

The  hydraulic  shaker  system  to  be  used  in  the  analysis  and  simula¬ 
tion  to  follow  consists  of  an  amplifier,  servovalve,  hydraulic  amplifier, 
piston  actuator,  and  feedback  transducers.  A  block  diagram  of  the  sys¬ 
tem  is  shown  in  Fig.  2. 


Fig.  2  Block  Diagram  of  Hydraulic  Shaker  System 
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The  servoamplifier  is  used  to  sum  the  input  and  feedback  signals 
and  to  provide  adjustable  gain  in  the  forward  loop.  The  gains  associ¬ 
ated  with  the  feedback  signals  are  also  adjustable.  The  amplifier  pro¬ 
duces  a  current  in  the  coil  of  the  electrohydraulic  servovalve  which,  in 
turn,  controls  a  flow  to  the  hydraulic  amplifer  or  power  spool.  The  - 
power  spool  then  controls  (with  a  large  flow  capability)  the  flow  to  the 
piston  actuator. 

In  order  to  have  fine  prefire  control  in  establishing  a  zero  position 
for  the  rocket  motor,  the  shaker  system  is  connected  as  a  positional 
servosystem.  A  signal  proportional  to  the  movement  of  the  piston  with 
respect  to  the  actuator  housing  is  fed  back  and  compared  to  the  input 
signal.  Because  the  actuator  is  basically  an  integral  device,  a  signal 
proportional  to  the  flow  to  the  actuator  (position  of  the  power  spool)  is 
fed  back  to  add  damping  to  the  system. 

2,3.2  Frequency  Response  Characteristics 


The  frequency  response  of  a  hydraulic  shaker  system  is  usually 
given  in  terms  of  relative  output  acceleration  versus  frequency.  Rela¬ 
tive  acceleration,  G,  is  the  ratio  of  actual  output  acceleration  to  gravi¬ 
tational  acceleration,  g.  Because  the  output  acceleration  is  not  only  a 
function  of  the  input  signal,  but  also  a  function  of  load,  the  frequency 
response  is  represented  by  a  family  of  curves.  Each  curve  is  the 
response  for  a  single  value  of  low  resistance  and  massive  load.  A 
typical  straight-line  response  is  shown  in  Fig.  3  (Ref.  7). 


Fig.  3  Typical  Hydraulic  Shaker  Frequency  Response 


Region  A  has  a  slope  of  12  db  per  octave  and  is  generally  called  the 
region  of  maximum  stroke.  Over  this  region  the  actuator  excursions 
are  limited  only  by  the  design  limit  of  the  actuator;  reserve  pressure  is 
available  for  accelerating  the  load,  and  the  servovalve  is  still  capable 
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of  supplying  more  flow  if  needed.  Region  B  has  a  slope  of  6  db  per 
octave  and  is  referred  to  as  the  valve  flow  limit  region.  In  this  region 
of  operation,  limiting  velocity  of  the  piston  occurs  as  a  result  of  the 
limiting  flow  through  the  servovalve.  Region  C  is  the  region  of  con¬ 
stant  acceleration.  Here,  not  only  limiting  flow  but  an  acceleration 
limit  is  present  because  all  available  pressure  is  being  used  to  acceler¬ 
ate  the  load.  Other  break  points,  farther  out  in  the  frequency  response, 
occur  as  a  result  of  compressibility  effects  and  servovalve  response. 
For  large  loads,  region  B  may  not  exist,  and  operation  may  only  in¬ 
volve  the  stroke  and  acceleration  limit  regions. 


For  a  fixed  load  condition  and  a  constant -amplitude  input  signal, 
the  shaker  system  can  be  represented  by  the  following  transfer  function: 


G(s) 


_ K_£ _ 

(T.s+1)  (Tbs+1>  (Tcs+  1) 


(1) 


A  straight-line  logarithmic  magnitude  plot  of  this  function  is  shown  in 
Fig.  4. 


lOg  CD  — 


Fig,  4  Response  for  o  Single  Constant  Load 


The  first  break  point,  wa,  is  fixed  by  the  servovalve  design.  In¬ 
vestigation  of  typical  commercial  response  curves  indicates  that  ua 
varies  inversely  with  load,  and  <jc  varies  inversely  with  the  load  raised 
to  some  power  in  fitting  the  transfer  function  to  the  frequency  response 
curves  for  other  loads. 
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From  the  standpoint  of  analysis  and  simulation  of  the  entire  thrust 
stand  system,  Eq.  (1)  is  of  little  value.  The  load  on  the  hydraulic 
actuator  varies  drastically  during  the  firing  and  tailoff  transients.  A 
Linear  analysis  would  have  little  meaning,  and  the  problem  must  be 
approached  through  a  detailed  analysis  of  individual  components. 


SECTION  III 

SYSTEM  ANALYSIS  AND  COMPONENT  SELECTION 


This  section  will  be  concerned  with  a  detailed  analysis  of  the  thrust 
system  to  be  simulated.  Describing  equations  will  be  developed  and 
specific  values  or  ranges  of  values  will  be  assigned  to  the  coefficients. 
In  the  case  of  commercial  components  such  as  in  the  hydraulic  servo- 
system,  specific  models  will  be  selected  and  the  manufacturers  litera¬ 
ture  used  as  a  guide  in  determining  transfer  functions. 


3.1  THRUST  BUTT 


With  reference  to  Fig.  1,  the  thrust  butt  is  to  be  represented  by  a 
linear  second-order  system.  The  force  equation  describing  the 
dynamics  of  this  system  is 

Mb^  Xb  +  Cb  sXb  +  KbXb  =  Fb  (2) 

where  is  the  force  transmitted  by  the  hydraulic  actuator.  The  mass, 
M5,  will  be  assumed  to  include  the  mass  of  the  actuator  housing, 
hydraulic  amplifier,  and  servovalve. 


The  undamped  natural  frequency  of  this  system  is 


and  the  damping  is 

C  _  I  cb 
2  yj  KbMb 

Under  static  load  conditions  the  steady-state  displacement  is 


(3) 


(4) 


(5) 
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Typical  values  of  40  cps  undamped  natural  frequency,  0.  008  damping 
ratio,  and  a  displacement  of  0.  04  in.  per  100,  000  lbf  leads  to  the 
following  assigned  values  for  the  coefficients  of  Eq.  (2): 

Kb  =  2.5  x  106  lb/in. 

Mb  =  39.61  lb  sec  /in. 

Cb  =  159.2  lb  sec/in. 


3.2  ENGINE  AND  LOAD.CELL  COMBINATION 

The  mass  of  the  rocket  engine  and  the  compliance  and  damping  of 
the  load  cell  form  a  mechanically  tuned  system.  The  force  equation 
describing  the  dynamics  of  this  system  is 

Mm  s’Xm  -  Cm  5  (XB1  -  X.)  +  Km  (Xm  -  Xa)  =  Tm  (6) 

The  force,  Tm,  is  the  thrust  produced  by  the  rocket  motor.  Depending 
on  the  size  of  the  engine  and  the  load  cell,  the  damped  natural  frequency 
(for  Xa  =  o)  generally  falls  in  the  range  from  30  to  50  cps.  The  fre¬ 
quency  may  vary  during  the  firing  of  solid-propellant  rocket  engines 
because  of  the  mass  change  caused  by  burning.  So  that  the  simulation 
will  represent  a  worst  case,  the  undamped  natural  frequency  will  pur¬ 
posely  be  set  to  44  cps  (slightly  higher  than  the  thrust  butt  frequency) 
and  the  damping  to  0.  008.  The  static  displacement  will  be  set  for  0.1  in. 
per  100,  000  lbf.  The  following  assigned  coefficients  will  produce  a  sys¬ 
tem  with  these  characteristics: 

Km  «  1  x  10*  lb/in. 

Mm  =  13.1  lb  secVin. 

Cm  =  57.92  lb  sec/in. 


The  force  as  seen  by  the  load  cell  is 

I'm  =  Cm  a(X[n  —  Xa)  +  Km  —  Xa) 
so  that  Eq,  (3)  may  be  written  as  follows: 

Mm  B2Xm  +  Fm  =  Tni 


(7) 

(8) 


3.3  HYDRAULIC  ACTUATOR 

The  hydraulic  actuator  of  Fig.  1  is  shown  in  more  detail  in  Fig.  5. 
For  this  analysis  it  will  be  assumed  that  the  active  areas  of  the  piston 
are  equal,  and  that  the  mass  of  the  piston,  Mp,  includes  that  of  the 
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piston  rod.  Also,  the  piston  seals  will  be  assumed  to  have  no  leakage. 
The  friction  will  be  of  the  viscous  type  with  no  stiction. 


Fig.  5  Simplified  Hydraulic  Piston  Actuator 


3.3,1  Development  of  Describing  Equations 

The  actuator  force  transmitted  to  the  thrust  butt  is 

Fb  -  PaAp  -  Pa'Ap  -  Cp  s(Xb  -  Xa)  (9) 

The  force  balance  equation  for  the  piston  is 

Fh  =  F™  -  Mp  ssXa  (10) 


The  pressures  on  either  side  of  the  piston  are  functions  of  the  posi¬ 
tion  of  the  piston  and  of  the  fluids  present.  Stated  mathematically  the 
pressures  are  defined  as  follows: 

Pa  =  <tf(Xb  -  Xa),  WJ  (11) 

Pa'  =  S((Xb  -  Xa),  W,1  (12) 

With  time  the  dependent  variable,  the  total  derivative  of  Pa  is 


dP. 

^(xb  -  X.") 


d(xt  -  X,)  <3P , 

dt  diV, 


di 

Ub~  X.) 


(13) 


From  the  definition  of  bulk  modulus,  jS,  partial  derivatives  of  Eq.  (13) 
may  be  expressed  as  follows: 


dPa  /3A„ 

d(xh  -  x.)  v 


(14) 
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and 


t?P.  _  JS 
dwx  *v 


(15) 


where  V  is  the  fluid  volume  in  one  end  with  the  piston  centered, and  w  is 
the  specific  weight  of  the  fluid.  Substituting  Eqs.  (14)  and  (15)  into 
Eq.  (13)  and  taking  the  Laplace  transform 

sPa  =  -iLsW,  -  ^Ps<Xb  -  Xa)  (16) 


In  like  manner  Eq.  (12)  may  be  written 


sP;  =  - 


wV 


sWa 


/3A»  tx 
—  sCXb  - 


Xa) 


(17) 


Equations  (9),  (10),  (16),  and  (17)  describe  the  dynamics  of  the  actuator 
in  terms  of  system  variables.  To  arrive  at  an  expression  that  will 
more  qualitatively  indicate  response,  the  undamped  natural  frequency 
and  damping  of  the  actuator  must  be  determined. 


Integrating  Eqs.  (16)  and  (17)  with  respect  to  time,  and  multi¬ 
plying  each  by  the  active  area  of  the  piston,  yields 

/3a> 


PaAp  = 


Pa'Ap  = 


j8ApWs 


wV 


(Xb  -  Xa> 


<xb  -  x.) 


(18) 

(19) 


Substituting  Eqs.  (18)  and  (19)  into  Eq,  (9) 
/3apw,  $Ap 


or 


/3ApW2  0A* 


(Xb  -  X.)  ♦  - 


=  Fm  -  MpSsXa  +  Cp5  (Xb  -  Xa) 


(Xb  -  xa) 


SJXa  + 


M 


■sXn  + 


p 

iSApW, 

MpwV 


2j3Aa„ 

VMp 

iSAp 


Xm  P  v 

g  =  ■  "I"  8A  V 


v  0APW,  jSApa  v 
Ab  ~  -  v»,  +  Hu  Ab 


VMp  D  wVMp  VMp 


(20) 


(21) 


From  Eq.  (21)  the  undamped  natural  frequency  of  the  actuator  is 


i 

f„  =  — 

2it  yj 


mi 

VM_ 


(22) 
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and  the  damping  coefficient  is  calculated  as  follows: 


2£«„  = 


cp 


M 


p 


c  = 


VMp 

2^Ap2 


(23) 


for  constant  Xfcj,  W^,  and  W2.  The  frequencies  observed  with  a  simple 
actuator  and  massive  load  cannot  be  entirely  predicted  with  Eq.  (22) 
for  the  following  reasons : 

(1) 


(2) 


(3) 


For  simulation  purposes  the  primary  use  of  Eq,  (23)  is  to  determine 
the  value  of  Cp  for  a  desired  value  of  damping.  This  will  be  needed 
when  numerical  values  are  assigned  to  the  coefficients  of  the  describing 
equations. 


The  analysis  leading  to  Eq.  (22)  did  not  account  for  stiction. 
Coulomb  friction,  and  leakage  past  the  piston.  Even  if  these 
nonlinearities  did  not  exist,  the  damped  natural  frequency 
would  be  lower  than  the  undamped  natural  frequency  given  by 
Eq.  (13)  as  a  function  of  damping. 

Under  some  conditions  of  excitation,  cavitation  occurs,  alter¬ 
nately  making  Pa  and  approximately  zero.  This  effect  tends 

to  reduce  the  natural  frequency  by 

Oscillations  about  an  off-center  point  where  the  two  fluid 
volumes  are  not  equal  causes  the  fluid  compliance  to  be  greater 
in  one  direction  than  in  the  other.  This  causes  the  observed 
frequency  to  be  higher  and  more  distorted. 


1 


3.3.2  Actuator  Specifications 

For  this  simulation  study  an  actuator  with  the  following  specifica¬ 
tions  was  selected  (Ref.  8): 


Maximum  dynamic  force 
Maximum  stroke 
Maximum  velocity 
Piston  weight 
Active  piston  area 
Total  fluid  volume 


100,  000  lb 
2  in.  DA 
9  in,  / sec 
220  lb 
41.67  in.2 
102.  0  in.3 
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The  hydraulic  fluid  will  be  assumed  to  have  a  bulk  modulus  of 
250,  000  psi  and  a  specific  weight  of  53  Ib/ft^.  Using  Eq.  (23),  the 
value  of  Cp  for  a  damping  coefficient  of  0.  2  and  other  coefficients  in 
agreement  with  the  above  specifications  is  as  follows: 


Cp  =  2Ap  C 


_  2(41*67)  (0  2)  /  x  106)  (220/386) 


102/2 


=  1245.0  lb  sec/in. 


3.4  SERVOVALVE  AND  HYDRAULIC  AMPLIFIER 


With  reference  to  Fig.  2,  the  transfer  function  for  the  servovalve 
and  hydraulic  amplifier  to  be  used  in  the  simulation  study  are  as 
follows  (Ref.  9): 


where 


ei>,  “  3768  radians/sec 
cuj  =  1884  radians/sec 
Ci  =  0«2o 


The  gains,  and  K2,  as  well  as  the  gains  of  the  servoamplifier  and 
feedback  transducers  will  not  be  specified  since  these  will  be  adjusted 
for  optimum  performance  prior  to  the  collection  of  simulation  data. 


The  flow  specification  for  the  selected  valve  is  100  gpm  at 
2500-psi  pressure  drop.  Using  the  general  flow  equation  for  flow 
through  a  restriction. 


or 


and 


max 
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In  terms  of  weight  flow,  the  flow  equation  can  be  written  as  follows: 


sW  =  K'XV  (0.00223)  w  -\/ap" 
Substituting  K'  from  Eq.  (27). 


(28) 


sW  =  0.00446  w 


./Xp 

x.l  V 


(29) 


^  v 

The  quantity,  — — —  ,  is  the  normalized  value  of  Xv  and  will  be  de- 

max 


fined  as  follows: 


£  ^vn,  (“1  £  ^vn  <  +  1) 


With  reference  to  Fig.  5,  the  flows 
using  Eq.  (29)  as  follows: 

SW,  =  0.2364  X 

sW3  =  0.2364  xvn 

and 


dW- 

dt 


and 


dW< 


dt 


will  be  defined 


vn  "\j  ^a  “  Pa 


%/p:  -  p»  \ 


for  Xvn  Positive 


sff,  =  0.2364  Xvn  yj  Pa  -  P0  J 

Vp*  - 


for  XVD  Negative 


sWa  =  0.2364  Xv 


(30) 

(31) 

(32) 

(33) 


The  simulation  circuit  will  be  designed  to  select  the  appropriate  flow 
equation  by  sensing  the  sign  of  Xvn. 


SECTION  IV 

SIMULATION  RESULTS  AND  DISCUSSION 


The  equations  resulting  from  the  analysis  of  the  thrust  stand  sys¬ 
tem  were  magnitude  and  time  scaled  (Appendix  I).  The  scaled  equations 
were  then  programed  on  an  analog  computer  to  form  an  analog  simula¬ 
tion  of  the  system.  Two  computer  circuits  were  used  for  the  study: 

(1)  simulation  of  the  basic  stand  without  the  servosystem,  and  (2)  sim¬ 
ulation  of  the  thrust  stand  system  with  the  servo  installed  (Appendix  II). 
Both  simulation  circuits  were  thoroughly  checked  for  program  and 
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wiring  errors.  They  were  then  subjected  to  a  series  of  operational 
checks  to  confirm  static  and  dynamic  performance. 

This  final  section  of  this  report  will  be  concerned  with  an  investi¬ 
gation  of  the  performance  of  the  thrust  stand  for  different  configurations. 
Both  the  time  and  frequency  domains  will  be  investigated  and  discussed 
in  the  material  that  follows . 


4.1  CONFIGURATION  1 


In  order  to  establish  a  reference  for  comparison  purposes  the  basic 
thrust  stand  without  the  servo  was  investigated  first.  Figure  6  is  a 
diagrammatic  representation  of  this  system. 


K 


b 


K 

m 


A  -  -  -  rr _ 

Mb 

_ UJi _ 

- rr _ _ 

0 

J  LL 

cb 

C 

m 

Xm 

T 

m 


Fig.  6  Basic  Thrust  Stand  without  the  Servo  Installed 

The  equations  describing  the  dynamics  of  this  system  are  as  follows: 


Mbs*Xb  +  CbsXb  +  KbXb  = 

Fm 

(34) 

[j-jS^Xin  +  Cms(Xm  —  Xfo)  + 

-  Xb)  =  Tm 

(35) 

fm  =  Cmsixm  -  Xfc)  +  Km(Xm 

Xl 

X 

i 

(36) 

Using  the  same  values  for  the  constants  in  these  equations  as  were  used 
in  the  analysis,  this  system  was  subjected  to  several  tests  as  explained 
below. 

A  limited  ramp  function  was  used  for  the  thrust  function.  Several 
rise  times  varying  from  0.  2  to  0.  01  sec  were  used  with  the  limit  level 
in  every  case  set  to  20,  000  lb.  The  applied  thrust  functions,  Tm,  and 
the  load-cell  responses,  Fm,  were  recorded  as  shown  in  Fig.  7.  The 
ramp  function  was  then  replaced  with  a  1,  000-lb  peak  sinusoidal  thrust 
function.  The  frequency  was  varied  to  obtain  the  load-cell  and  motor 
position  frequency  responses  shown  in  Fig.  8. 
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In  Fig.  8  the  high  Q  resonance  at  30  cps  checks  very  closely  with 
the  damped  natural  frequency  in  Fig.  7  for  the  fastest  rise  times.  The 
low  frequency  shift  in  the  level  of  Fm  in  Fig.  7(e)  occurring  at  approxi¬ 
mately  5.  5  cps  is  confirmed  in  Appendix  III. 

In  order  to  check  the  resonant  frequencies  appearing  in  Fig.  8, 

Eqs.  (34),  (35),  and  ( 36)  were  solved  simultaneously  letting 

Cm  =  Cb  =  0  to  yield  the  following  auxiliary  equation  (set  Fm  «=  0  and 

substitute  r  for  s) 


<  ^  KmKb  +  KbMm  +  Kn,Mnl  ,  KmKb  _  Q  (37) 

Since  the  complex  roots  of  a  polynomial  with  real  coefficients  occur  in 
conjugate  pairs  and  for  this  check  the  real  parts  are  known  to  be  zero 
(Cm  =  Cb  =  0),  the  roots  of  Eq.  (37)  provide  only  two  frequencies.  Sub¬ 
stituting  the  assigned  numerical  values  for  the  constants  in  Eq.  (37) 
and  factoring,  the  roots  are  found  to  be  31  and  56.  7.  These  values 
check  very  closely  to  the  resonant  frequencies  of  30  and  55  cps  shown 
in  Fig.  8.  The  resonant  frequencies  should  be  slightly  below  the  un¬ 
damped  natural  frequencies  calculated  above. 


4.2  CONFIGURATION  2 

The  thrust  stand  system  with  the  servo  installed  is  considered  in 
this  investigation.  The  servo  is  connected  as  a  positional  servosystem 
where  the  command  signal  is  furnished  by  the  manual  set-point  poten¬ 
tiometer,  and  the  piston  position,  Xb  -  Xa,  is  the  feedback  signal. 
Because  the  piston  is  basically  an  integral  device  (for  light  loads  and 
negligible  compression  of  the  hydraulic  fluid),  the  position  of  the  power 
spool  is  proportional  to  the  piston  velocity. '  A  proportional  amount  of 
the  power  spool  position,  Xv,  is  fed  back  to  provide  damping  to  the 
servosystem. 

Starting  with  a  very  small  damping  signal,  the  servo  gain  was  ad¬ 
justed  for  a  stable  but  slightly  oscillatory  response  for  a  small  step 
change  in  the  manual  set-point  signal  (10  percent).  The  responses  of 
-  Xa,  Xv,  and  Fm  were  observed  on  the  repetitive  display  scope  of 
the  computer.  With  each  increase  in  the  servo  gain,  the  damping  was 
increased,  accordingly.  A  point  was  finally  reached  where  a  more 
stable  operation  could  not  be  gotten  by  increasing  the  damping  adjust¬ 
ment.  At  this  point  the  servo  gain  was  lowered  slightly  and  the  damp¬ 
ing  adjusted  for  optimum  response.  The  data  of  Fig.  9  were  then 
recorded. 
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The  highly  oscillatory  response  of  Xv  was  investigated  and  found  to 
be  close  to  the  damped  natural  frequency  of  the  servovalve.  The  damp¬ 
ing  on  the  valve  was  increased  from  0.  25  to  0. 7  and  the  data  of  Fig.  10 
were  recorded  to  show  the  effect.  Comparing  Figs.  9  and  10,  the 
response  of  Xv  was  greatly  affected,  but  that  of  Xb  -  Xa  and  Fm  were 
affected  only  slightly.  Because  the  damping  was  originally  specified  to 
be  0.  25,  and  this  would  not  normally  be  an  adjustment  on  the  servovalve, 
this  value  was  used  for  the  rest  of  the  study. 


-5,000  I - 1 - 1 _ I - 1 _ I _ 1  1  I  i  I 

0  0.5  1.0 


Time,  sec 

Fig.  10  Response  in  Manual  Mode  for  Step  Change  in  Servo  Command  Signal  with 
Pilot  Valve  Damping  Set  tor  £  -  0.7 
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In  order  to  compare  the  performance  of  this  configuration  with  that 
of  Configuration  1,  the  system  was  tested  with  the  same  thrust  functions 
(ramp  and  sinusoidal)  to  obtain  the  data  shown  in  Figs.  11  and  12.  Com¬ 
paring  Figs.  7  and  11,  the  load  cell  responses  are  almost  identical, 
except  that  in  Fig.  7  the  damping  is  slightly  less.  The  frequency 
responses  of  Figs.  8  and  12  are  also  very  much  alike  except  for  the 
higher  Q  resonances  in  Fig.  8;  this  correlates  well  with  the  time 
responses  since  lower  damping  and  higher  Q  are  directly  related. 


4.3  CONFIGURATION  3 


This  investigation  will  be  concerned  with  the  automatic  mode  of 
operation.  In  this  mode  the  command  signal  to  the  servo  (replacing  the 
manual  set-point  signal)  is  derived  from  live  measurements  from  the 
rocket  motor  and  thrust  stand.  For  reasons  explained  below,  the  com¬ 
posite  signals  considered  practical  for  this  purpose  are  the  load-cell 
output,  Fm,  and  the  rocket  motor  velocity  and  acceleration,  sXm  and 
s^Xm,  respectively.  The  gain  and  damping  determined  in  Configura¬ 
tion  2  for  the  servosystem  is  also  used  here  since  the  servosystem  is 
still  used  to  position  the  piston  actuator. 

The  composition  of  the  automatic  feedback  signal  was  determined 
primarily  on  the  basis  of  the  practical  aspects  of  the  problem.  Gener¬ 
ally,  Fm,  sXm,  and  s2Xm  are  readily  available  from  single  instru¬ 
ments  with  only  minor  instrumentation  problems.  Direct  measurement 
of  motor  position  is  practically  impossible  in  a  test  cell  because  of  the 
absence  of  a  fixed  reference.  Measurement  of  motor  position  by  inte¬ 
gration  of  the  velocity  is  accompanied  by  noise  and  drift  problems. 


In  order  to  satisfy  the  performance  goal  of  zero  motor  position 
during  the  steady  portion  of  the  firing,  a  signal  proportional  to  the 
steady  level  of  thrust  was  needed.  This  was  provided  by  the  load  cell 
output  since  Tm  is  equal  to  Fm  in  steady  state.  For  zero  rocket  motor 
position  in  steady  state,  the  position  of  the  piston  is  related  to  the  load¬ 
cell  output  by  the  following  equation: 


Xb 


Xa  =  f-i—  +  -L\  Fj, 

K  J 

-  KtFm 


(38) 
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Fig.  12  Frequency  Response  of  Thrust  Stand  System  in  the  Manual  Mode 
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Feeding  KtFm  to  the  servoamplifier  in  the  simulation  caused  the 
system  to  become  unstable.  This  condition  was  brought  under  control 
by  introducing  a  proportional  amount  of  motor  acceleration.  Velocity 
was  also  added  with  some  improvement.  After  optimizing  the  gains 
associated  with  velocity  and  acceleration  in  the  feedback  path,  the  gain 
on  the  velocity  term  was  so  small  that  this  signal  was  removed  before 
any  data  were  recorded. 

Figures  13  and  14  show  the  response  of  the  thrust  system. 

Figures  13a  and  14a  show  the  response  for  ramp  thrust  functions  of 
0.  2-  and  0.  01-sec  rise  times  in  the  manual  mode  with  the  manual  set- 
point  potentiometer  set  to  zero.  Figures  13b  and  14b  show  the  responses 
to  the  same  thrust  functions  in  the  automatic  mode  with  the  active  feed¬ 
back  set  to  KtFm  plus  an  optimized  amount  of  s2xm.  The  limit  level  of 
the  thrust  function  in  every  case  was  20,  000  lb.  With  the  same  feedback, 
the  ramp  thrust  function  was  replaced  with  a  sinusoidal  thrust  function 
(1000-lb  peak)  and  the  frequency  response  data  of  Fig.  15  were  collected. 
The  following  characteristics  and  improvements  may  be  seen  from  these 
data  and  comparisons  to  previously  discussed  data: 

1.  A  high  degree  of  damping  has  been  introduced  into  the  system 
in  the  automatic  mode  as  seen  by  the  time  responses  of 
Figs.  13  and  14.  In  the  frequency  domain  this  damping  has 
greatly  reduced  the  Q  of  the  resonances.  Notice,  too,  that  the 
resonant  frequencies  have  been  slightly  shifted. 

2.  In  the  low  frequency  region  of  both  Figs.  8  and  12  the  response 
of  the  motor  position  is  flat,  indicating  that  over  this  frequency 
range  the  position  is  directly  proportional  to  the  thrust  ampli¬ 
tude.  In  Fig.  15  the  closed  loop  control  has  provided  a  high 
degree  of  attenuation,  especially  in  the  low  frequency  region  of 
Xm.  The  motor  position  is  not  only  proportional  to  thrust  am¬ 
plitude  but  also  proportional  to  frequency;  the  slope  is  approxi¬ 
mately  20  db/decade.  The  displacement  at  zero  frequency 
(steady-state)  is  zero  as  shown  in  Figs.  13b  and  14b. 

3.  Comparing  the  frequency  responses  of  the  basic  stand  (Fig.  8) 
with  those  of  Fig.  15,  the  automatic  mode  has  extended  the 
frequency  response  of  the  load  cell  and  reduced  that  of  the 
rocket  motor  position.  Poor  response  in  the  rocket  motor 
position  at  higher  frequencies  is  favorable  toward  reducing 
the  offset  during  the  firing  transient. 
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4.  In  Fig.  8,  the  resonance  at  235  cps  was  investigated  and  found 
to  be  caused  by  the  resonant  frequency  of  the  servomotor. 

This  is  slightly  less  than  the  damped  natural  frequency  of 

260  cps  observed  in  the  response  of  Xv  in  Fig.  9.  The  damp¬ 
ing  of  the  servomotor  was  increased  from  0.  25  to  0.  5,  and  the 
frequency  response  of  Fm  checked  again.  The  response  below 
100  cps  remained  unchanged  but  that  in  the  high  frequency 
range  followed  the  dotted  portion  of  the  curve  shown  in  Fig.  15. 

5.  For  rocket  motors  where  there  is  no  steady  portion  of  the 
firing  but  a  steady  drop  in  thrust,  the  position  of  the  rocket 
motor  will  take  on  a  value  proportional  to  the  slope  of  the 


27 


AEDC-TR-67-232 


thrust  (lb /sec).  This  is  obvious  from  the  time  response  of 
Xm  in  Fig.  13b  during  the  rise  time  of  the  thrust  function. 

This  could  be  corrected  by  introducing  a  small  amount  of  the 
integral  of  position  as  part  of  the  feedback  signal.  Another 
way  might  be  to  use  a  signal  which  is  composed  of  propor¬ 
tional  plus  derivative  of  the  load-cell  output.  In  any  event, 
the  correcting  signal  would  be  small  (dependent  on  the  slope 
of  the  thrust  function)  in  most  cases,  which  would  favor  the 
problem  of  working  with  integrators  or  adding  lead  networks. 

The  ramp  and  sinusoidal  thrust  functions  used  thus  far  have  been 
small  enough  to  ensure  operation  within  limiting  velocity  limits.  As 
the  thrust  level  goes  up,  however,  the  piston  is  required  to  move 
faster  and  over  a  greater  amplitude,  and  limiting  velocity  begins  to 
affect  the  response.  Figure  16  shows  the  response  of  motor  position 
and  the  load  cell  output  for  various  levels  of  applied  thrust.  The  rise 
time  of  the  thrust  functions  is  0.  01  sec  in  every  case.  Remembering 
that  the  normalized  value  of  Xv  is  Xvn,  which  ranges  in  value  from  -1 
to  1,  the  exact  point  at  which  limiting  velocity  occurs  can  be  observed. 

For  thrust  levels  up  to  60,  000  lb  the  responses  of  Xm  and  Fm 
appear  nearly  linear.  Limiting  velocity  is  just  reached  when  the  thrust 
function  is  increased  to  70,  000  lb.  Beyond  this  value  the  flow  limit  is 
reached  more  than  once,  and  the  response  of  Fm  becomes  noticeably 
affected.  For  the  80,  000-lb  thrust  function  the  pressure  P^  dropped  to 
9.  5  psia.  Then,  for  90,  000  lb  thrust,  dropped  to  4  psia.  This  is 
the  region  at  which  cavitation  occurs  and  the  equations  concerning  flow 
become  invalid. 


SECTION  V 
CONCLUSIONS 


The  results  of  this  simulation  study  indicate  that  the  dynamic 
response  of  existing  thrust  stands  can  be  greatly  improved  using  a  high 
quality  servosystem.  A  system  similar  to  Configuration  3  would  pro¬ 
vide  the  following  improvements  in  comparison  to  a  basic  stand. 

1.  A  high  degree  of  system  damping  would  be  introduced. 

2.  Immediately  after  the  firing  transient  the  rocket  motor  would 
return  to  its  initial  prefire  position.  The  motor  would  remain 
at  this  position  even  for  moderately  rapid  changes  in  thrust 
level  because  of  the  fast  tracking  capability  of  the  servosystem. 
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Xm,  In. 


20,.000-lb  Applied  Thrust 


ai 

0 


b.  60,000-lb  Applied  Thrust 


e.  70,000-lb  Applied  Thrust 


100,000 
50,000 

ol - 1 - 1 

0  0.2 


o  0.2 


Rise  Time,  sec 


e.  90,000-lb  Applied  Thrust 


Fig.  16  Response  of  Thrust  Stand  in  the  Automatic  Mode  to  Varying  Thrust  Levels 
with  0.01. sec  Rise  Times 
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3.  Optimization  of  the  control  circuitry  would  greatly  decrease 
the  offset  of  the  motor  position  during  the  firing  transient. 

4.  Concurrent  with  the  improvements  described  above,  the 
acceleration  of  the  rocket  would  be  greatly  decreased  during 
the  firing  and  tailoff  transients. 
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APPENDIX  I 

TIME  AND  AMPLITUDE  SCALING  OF  SYSTEM  EQUATIONS 


The  equations  resulting  from  the  analysis  of  the  thrust  system  are 
summarized  below.  Some  of  the  equations  from  the  text  have  been 
arranged  below  so  that  they  express  the  quantity  to  be  solved  for  by  the 
computer. 


t'b  _ v.  Kb 


sJXb 

F'm  =  CmfitXfn  —  Xa)  +  Km  (Xm  *"  Xa) 

Fb 

s2xa  =  (rm  -  Fb) 


«b  ~  MrsXb  -  nr Xb 


J_(T„  -  Fm) 

it  v  1  m  1  nv 


=  (Pa  -  Pi)  Ap  -  Cph  (Xb  -  Xa) 


sPa  - 


^_sWl  -  ^s(Xb  -  xa> 


aP«  “  "  w  sW»  +  HT  s(Xb  -  X**> 


Xp(±-  +  l)  =  , 

p  \o)j  )  sJ  -  2  4jW3  b  i-  Wj 


x  _  JL  x 

Ay  ~  Ap 


\ 

SW,  =  0.00446  wXvn  J  Ps  -  Pa  / 
■W,  =  0.00446  wXvn  yj  P;  -  P0^ 


’>  for  Xvn  Positive 


sWj  =  0.00446  wXvn  yj  PB  -  Pa 
sW,  -  0.00446  wXvn  yj  Ps  -  P'  \ 


1  for  X  vn  Negative 


(2) 

(7) 

(8) 

(9) 

(10) 

(16) 

(17) 

(24) 

(25) 

(30) 

(31) 

(32) 

(33) 
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The  coefficients  which  apply  to  the  above  equations  are  summarized 
below: 


K i,  =  2.5  x  104  lb/in. 
M(,  =  39.61  Hi  /in, 
C |,  =  159.2  lb  Hpc/in. 
km  =  1  x  10*  lb/in. 

Mm  =  13.1  lb  secVin. 
Cm  -  57.92  lb  sec/in. 
Ap  =  41.67  in.J 
V  =  51.0  in.1 


w,  =  3768.0  rad  ians/sec 
oi,  =  1884.0  radians/sec 
£  =  0.25 
j8  =  0,25  x  10*  psi 
Mp  =  0.5699  lb  secVin. 
Cp  =  1245.0  lb  sec/in. 
w  =  53.0  lb/ft’ 


Substituting  the  above  coefficients  into  the  system  equations  and 

scaling  the  system  variables  so  that  the  quantities,  — ,  —  , 

L io3J  L 103J  L30j 

[lOO  w]  ,  and  [lOO  Xyn]  appear,  the  equations  become 


0.0524  X  10*  4.02sjio1Xbj  -  0.06311  x 


57.92  x  10*s|^103Xin  -  WX^J-  1.0  |lOsXm  - 


0.07633  x  10*  pJS-  _  Isl 

LlO1  10*J 


1.25  ?-±-~  -  0.001245s  jlO’Xb  -  WX.J 


1.755  x  10* 


53.27s  [l00W^  -  6.810s  [l0}Xb  -  l^Xa] 


-53.27  s  JlOOwjJ  -  6.810  sjlO’Xb  -  103XaJ 


(2) 

(7) 

(8) 


(9) 


(10) 

(16) 


(17) 
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Xp(2.654  xir8  +  l)  = 


In  order  to  facilitate  the  recording  of  simulation  data  on  an  X-Y 
recorder,  the  above  scaled  system  equations  are  slowed  down  in  time 
by  a  factor  of  100.  This  is  accomplished  by  letting  s  in  the  above  equa¬ 
tions  be  100s.  In  the  time  domain  this  is  equivalent  to  letting 

-f  =  100  where  r  =  100t 
dt  dr 


The  thrust  function  Tm  which  is  a  limited  ramp  function  with  rise 
times  ranging  from  0.  01  to  0.2  sec  will  be  time  scaled  accordingly  so 
that  the  rise  times  from  1.  0  to  20  sec  will  be  used. 

The  final  amplitude  and  time  scaled  system  equations  are  as 
follows : 


s3  |l0JXb]  =  2-524  [Tli]  -0.0402e  ^03Xb]  -  6.311 

fe]  =  °-005792  s  [1(PXn.  “  10>Xa]  +  1*°  [105Xn,  “ 

s*  [lO=X^]  =  7.633  [Tm  -J™  ] 

EtS>]  -  L2S  .  [lO>Xb  -  10>X.] 


(2) 

(7) 

(8) 

(9) 


35 


AEDC-TR-67-232 


s’  [lO’Xa]  =  175.5  [Fm 

a  [■?§]  =  53-27  s  (j00ffJ]  -  e.siOs  [lO’Xb  -  10^] 
SB0  =  -53,27s^OOW]  +  6.810  sjlO’Xb  -  10’XaJ 


X  p(0.02654  s  +  1) 


KE _ 

s’  +  9.425  +  354.9 


=  0.01295 


=  0.01295 


=  0.01295 


=  0.01295 


(10) 

(16) 

(17) 

(24) 

(25) 

(30) 

(31) 

(32) 

(33) 
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APPENDIX  II 
SIMULATION  CIRCUITS 


The  simulation  circuit  used  for  studying  Configuration  1  is  shown 
in  Fig.  II- 1.  Equations  (2),  (7),  and  (8)  were  used  for  this  study  by  letting 
Xa  =  and  F^  =  Fm.  The  scaled  versions  of  these  equations  are  shown 
at  the  end  of  Appendix  I. 


Pig.  1 1-1  Simulation  of  Basic  Thrust  Stand 


Figure  II - 2  is  the  analog  circuit  used  for  studying  the  second  two 
configurations.  This  circuit  is  the  analog  program  for  Eqs.  (2) 
through  (33)  at  the  end  of  Appendix  I. 
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Fig.  M-2  Simulation  Circuit  of  Thrust  Stand  and  Hydraulic  Servosystem 
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APPENDIX  III 

LOW  FREQUENCY  SHIFT  IN  THE  RESPONSE  OF  Fm 


In  Fig.  7e  two  frequencies  are  noticeable  in  the  response  of  Fm; 
these  frequencies  are  approximately  30  and  5.5  cps.  The  30-cps  fre¬ 
quency  is  the  damped  natural  frequency  corresponding  to  the  calculated 
undamped  natural  frequency  of  31  cps.  The  5.5-cps  frequency  is  not 
easily  related  to  either  the  31-  or  56.  7-cps  undamped  natural  frequencies 
or  to  combinations  of  the  two.  Since  this  low  frequency  is  also  present 
in  the  response  of  the  rocket  motor  position  (Fig.  14e),  its  presence  in 
the  response  of  Fm  can  best  be  explained  by  directing  our  attention  to 
the  analytic  solution  for  Xm  for  the  undamped  system. 

Letting  Cb  =  Cm  =  0  in  Eqs.  (34),  (35),  and  (36),  and  solving  these 
equations  simultaneously  for  Xjn,  an  equation  results  having  a  general 
solution  of  the  form  (Ref,  10) 

X[u  =  A,  +  A,  Ai  Cos  A,  Cos  cjjt  (JU“  1) 

At  t  =  0,  Xm  =  0  and  at  periodic  time  intervals  later  the  cosine  terms 
alternately  add  and  subtract  to  produce  variation  in  amplitude  at  the  rate 
of  5.  5  cps.  This  is  confirmed  in  the  computer  investigation  that  follows. 

Equations  (34),  (35),  and  (36)  were  scaled  and  programed  on  the 
analog  computer  (Fig.  Ill  - 1)  for  the  case  of  a  step  change  in  Tm  from 
zero  to  20,  000  lb.  The  system  damping  was  also  set  to  zero.  Two  un¬ 
damped  second-order  systems  were  set  up  to  produce  10  cos  uit  and 
10  Cos  U2t  where  uj  =  2?r(3l)  and  u>2  =  2ir(56.  7).  The  magnitude  of 
Ai  +  A 2  was  known  since  this  corresponds  to  the  steady-state  value  of 
Xm  (0.  028  in. )  when  the  steady  thrust  is  20,  000  lb.  The  sum  of 
Aj  +  A2  -  Ai  Cos  wit  -  A2  Cos  U2*  (output  of  amplifier  24  in  Fig.  III-1-b) 
was  observed  on  the  repetitive  operation  display  scope  simultaneously 
with  Xm  from  the  stand  simulation  (output  of  amplifier  12,  Fig.  III-1-a). 
The  coefficients  Hi  and  K2  were  manually  varied  until  the  two  solutions 
coincided.  The  values  for  Ai  and  A2  were  determined  to  be  22.  75  and 
5.  25,  respectively.  Figure  III -2  shows  the  responses  of  Xm  and  the 
general  solution  for  these  values  of  Ai  and  A2.  The  cosine  functions 
are  plotted  to  show  their  values  in  relation  to  the  response  functions. 

By  experimenting  with  other  frequency  combinations  it  was  deter¬ 
mined  that  the  low  frequency  shift  can  be  predicted  with  the  following 
equation: 

t  =■  |af,  -  U  I 

where  a  =  fj/h  rounded  to  the  nearest  integer. 

In  the  above  case  a  =  56.  7/31  =  1.  83  and  rounded  becomes  2.  0.  The 
low  frequency  calculated  using  this  criterion  is  5.  3  cps.  This  checks 
well  with  the  observed  frequency  of  5.  5  cps. 
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b.  Circuit  for  Generating  General  Solution 
Fig.  ill-1  Low  Frequency  Component  Investigation 
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